On the generalized Landau-Zener problem. 
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During the adiabatic time evolution levels crossing violates the adiabaticity and makes transitions 
between levels possible. Conventionally only two energy levels cross simultaneously. The transition 
probabilities for this case were found by Landau and Zener [jjj , ||. However, the multilevel crossing 
happens systematically rather than occasionally if the Hamiltonian possesses a special symmetry. 
The simplest physical realization of the multilevel crossing are the Zeeman multiplet in a varying 
magnetic field and an electron in one-dimensional chain driven by the time-dependent electric field. 
We present asymptotics of the transition amplitudes for these kinds of the multilevel crossing. They 
are based on an exact solution for a model n-state, time-dependent Schrodinger equation. 



Landau-Zener theory is one of the most important and 
influential results in non-stationary quantum mechanics. 
It has numerous applications ranging from the molecu- 
lar dissociation ||, slow atomic and molecular collisions 
H to electron transfer in biomolecules ||. Recently it 
was successfully employed for the experimental determi- 
nation of very small tunnel splittings || and quantum 
phase interference |7j in molecular nanomagnets. 

Landau-Zener (LZ) theory pj|, Q deals with an adi- 
abatic process in a quantum system with discrete spec- 
trum under external bias. The adiabaticity is violated 
when a pair of levels move towards each other strongly 
enhancing the transitions between the two states. LZ 
determined the transition probabilities for the two-level 
crossing. 

LZ theory considers the crossing of more than two lev- 
els at one moment of time to be an unlikely process. 
However, in some systems such a crossing may occur sys- 
tematically stemming from high symmetry of the under- 
lying Hamiltonian. A realistic example is a multiplet of 
atomic electronic states with the total spin S or total ro- 
tational moment J larger than 1/2 in varying external 
magnetic or electric field. The Zeeman splitting between 
25 + 1 or 2 J + 1 levels regularly vanishes at nodes of the 
magnetic field. 

Another example is provided by localized electronic 
states in a symmetric crystal environment. This degen- 
eracy can be lifted by driving external fields, as in the 
previous example and by the Jahn- Teller effect. A sim- 
ilar degeneracy is displayed by a large spin placed in a 
crystal environment of high symmetry JTo| ]. For some 
specific Hamiltonians employed in nanomagnets theory 
the existence of multilevel crossing at specific values of 
parameters (components of magnetic field) was shown 
P|, ||. Though this degeneracy probably is unstable 
with respect to perturbations of the Hamiltonian, these 
perturbations are small for real nanomagnets as it was 
demonstrated in the experiment by Wernsdorfer et al. 



In this work we find analytically the transition ampli- 
tudes for special cases of multilevel crossing. For this 
special symmetry the transition amplitudes can be ex- 
pressed in terms of a two-level problem. Physically this 
case corresponds to transitions in the Zeeman multiplet 
placed into a strong, time-dependent magnetic field di- 
rected along a constant axis (z) and a weak, slowly vary- 
ing or constant magnetic field in the perpendicular direc- 
tion (x). 

The time-dependent Shrodinger equations for two lev- 
els are: 



ibi = Ei(t)bi + A(t)b 2 
ib 2 = E 2 (t)b 2 + A(t)6i 



(1) 



Near a crossing point the dependence of energy levels 
on time is approximately linear E a {t) = E a t; a = 1,2, 
whereas the non-diagonal matrix elements of the Hamil- 
tonian can be taken constants. In terms of new ampli- 
tudes ai.2 = e~ l ( E i+ E 2) t / 4 6 12 , after a time rescaling, 
equations (0) can be simplified as follows fl: 



idi = ta\J1 + 7122 
ia 2 = —ta 2 /2 + 701 



(2) 



where 7 = A/yCt, (l = E\ — E 2 . Eliminating a 2 
from these equations, we find the parabolic cylinder 
equation for ai(t). Its solution which has asymptotics 

ai ~ exp(— i| i7 2 ln|i|) (and a 2 = 0) at t — ► —00 

is the Weber function D_j 7 2 (e l7r ^t) whose asymptotics 
are well known pj| . The scattering matrix for the two- 
level system is conveniently written in terms of modi- 
fied amplitudes c\ — a\ exp(if); c 2 — a 2 exp(—if) where 



/ = V + t 2 In 1*1- It reads: 



exp(— 7T7 2 ) 
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Let us consider a system with the total spin 5 > 1/2 in 
constant magnetic field H x along x-direction and vary- 
ing with time, much larger in average field H z (t) along 
z-direction. Its time evolution is regulated by the Hamil- 
tonian: 



Hq — —h x S x 



h z {t)S z 



(4) 



where h a = g^,BH a ; a = x,z and S Z ,S X are the spin 
operators. In he vicinity of its node h z (t) can be approx- 
imated by a linear function h z (t) = h z t where h z is the 
time derivative of h z (t) taken at the node. After a proper 
rescaling of time and energy the Hamiltonian (|]) takes a 
following form: 



H = 2^S X + tS z 



(5) 



It depends on one dimensionless parameter 7 = ,f (the 

LZ parameter). This Hamilton operator belongs to the 
50(3) algebra. This fact allows to derive the evolution of 
an arbitrary spin 5 in a varying magnetic field from the 
solution of Schrodinger equation for spin 1/2 in the same 
field. The corresponding evolution operator Us(t) is an 
operator of rotation belonging to the group 5*0(3) and 
acting in its irreducible representation. Since the compo- 
sition law does not depend on a specific representation, 
the resulting evolution operator at a fixed moment of 
time represents the same rotation for any spin. Thus, 



the problem is reduced to expression of the rotation op- 
erator for spin 5 if it is known for spin 1/2. Note that 
the Hamiltonians (||J|) are most general for the case of 
2-level crossing. 

The multi-spinor technique is most appropriate for this 
purpose (see [||, ch. VIII). In general the spin S state can 
be represented as a direct symmetric product of 2S spin 
1/2 states: 



(S + m)l(S-my 
l S ' m > = V (25)! (l " 



■) (6) 



where each ket contains 5 + m spins up and S — m spins 
down and all permutations are performed. Let the SU (2) 
matrix rotating spin 1/2 states is: 



a b 

—b* a* 



(7) 



with I a. 1 2 + \b\ 2 = 1. Equivalently an individual spinor is 
transformed according to: 



|+) _o|+>-&*|-): 
I-) -&|+)+ a *|_) 



(8) 



The transformation for the state (g) can be obtained as 
the direct product of transformations (ph. 



V2S(2S-1)\ aS+m( _ 6 , )S _ m _ v 



{S-m- l)!(5 + m)! 



I S, m) 



(S + m)\(S — m )j_ a S+m^_^*\ ) S-m |g g\ _|_ 



(25)! 



25(25 - 1) 



S-\-m — l { 7* \S—rn 



(S + m-iy.(S-m)\ 



b *y-m b \s,S-l) 



A general matrix element of the rotation operator (m | 

Us I m') for the spin S representation is expressed in terms of a, b, a*, b* in the following way Jl2| |, (L 



(m I Us I m) 



(S + m')\(S - m')\ 
(5 + m)!(5-m)! 



1/2 



a m ' +m b m '- m P™^ m ' +m {2\a\ 2 - 1) 



(9) 



where P^ b (x) are the Jacobi polynomials [|T3|. The ma- 
trix elements possess the following symmetry proper- 
ties: (-m I Us I -m') = (-1)M+I m 'l(m | U s \ m')*, 
|(m I U s I m')\ = |(m' | U s \ m)\ = \{-m' \ U s \ -m)|. 

Equation (^) displays oscillations of matrix elements 
associated with oscillatory behavior of the Jacobi poly- 
nomials. The number of nodes N{m, m') of the matrix 
elements (m | Us \ rn!) can be determined geometrically 
as the number of the square shell to which it belongs in 



the square matrix. We accept the number of the external 
square with max(\m\, \ m'\) — 5 for zero and number is 
increasing when the square shell size is decreasing. An- 
alytically N(m,m!) = S — max(\m\, \m'\). Due to sym- 
metry several matrix elements (2, 4 or 8) become zero at 
the same value of \a\ or 7. Central matrix elements have 
maximal number of nodes (5 for integer spins, 5—1/2 
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for half- integer spins). 

The scattering matrix for our problem can be found 
from the general expression (O) by substitution: 



H„ 



nh(t)S nm + jt 



(14) 



a = cxp(— 7T7 2 ), b = — 



2^exp(^_ 



7 r(-i7 2 ) 



(10) 



which follows from comparison of equations (§J7|). To 
make our results more visual, we present explicitly the 
scattering matrices for spins 5=1 and 5 — 3/2: 
5=1. 



Ui -- 
S = 3/2. 



a 2 y/2ab b 2 

-V2ab* 2\a\ 2 - 1 y/2a*b 
b* 2 -V2a*b* a* 2 



(11) 



where t nm = 1 if |n — m\ = 1 and zero otherwise. The 
physical implementation of this model is an electron tun- 
neling in a discrete equidistant chain placed into a vary- 
ing homogeneous electric field directed along the chain. 
We have introduced a natural set of states \n > located 
at the n-th site. The tunneling between the sites is sup- 
pressed by the external field since it lifts the degenera- 
tion of energy. However, when the field becomes zero the 
tunneling rate grows. Near this point we accept the time 
dependence of the electric field as linear. Then, after 
rescaling of time the initial Hamiltonian (|l4| ) reads: 



+ g{3m,n+l + 5m,n-l) 



(15) 



(a 3 



3/2 



V3a 2 b 



VEab 2 



~V3a 2 b* (3|a| 2 -2)a (3|a| 2 - 1)6 V3a*b : 
V3ab* 2 -(3|a| 2 -l)6* (3|a| 2 - 2)a* V3a* 2 b 
\ -b* 3 V3a*b* 2 -V3a* 2 fe* a* 3 



\ 



where g = j/\/eEA. It corresponds to a system of 
Shrodinger equations: 



nta n + g(a„_i + a n +x) 



(16) 



(12) 

To save space we present only several transition probabil- 
ities for 5 = 2, all the rest can be found from symmetry 
properties (see above): W 2 ,2-k = k \(t-k)\ (l Q | 2 ) 4 ~ fc (l _ 
|a| 2 ) fe , k = 0...4; Wx,x = |a| 4 (4|a| 2 - 3) 2 ; W^i.o = 6|a| 2 (l - 

| 2 ) 2 (2|a| 2 - l) 2 ; Wi,- X = (1 - M 2 ) 2 (4|a| 2 - l) 2 ; W ,o = 
(6|a| 4 -6|a| 2 + l) 2 . 

Additionally we present transition probabilities for 
5=1: 



/ The problem is to find the scattering matrix for this sys- 
tem. In other words, we are looking for an asymptotic 
at t — > +oo of a solution a n {t) which obeys the initial 
condition |et„(i)| 2 = 5 n .n' at t — > — oo. 

Let introduce an auxiliary function u(ip, t) = 
S^L-oo CL n (t)e mv . System ( |ili| ) is equivalent to a fol- 
lowing equation in partial derivatives for u(ip,t): 



du ,du 

— +t— 1- 2ig cos ipu = 

at dip 



(17) 



W hl = < 
Wx -i = ( 1 - e 



Wx,o = 2(e 



-47T7' 



); 



"* 2 ; W 00 = fl-2e- 27r72W 



(13) 



Other probabilities can be readily found from symmetry 
relations. As a consequence of Jacobi polynomials oscil- 
lation, the matrix elements of inner squares have nodes 
at some special values of 7. Thus, clement (Ux)oa is zero 

at \a\ 2 = 1/2 or at 7 = ^ ± In 2 sa 0.332. Elements of 

the matrix (1/3/2)1/2, 1/2 = (P3/2)-x/2-x/2 become zero 



We should find a solution of this equation which obeys 
the initial condition: u{ip,t) — > exp[m'(— 4j- + (p)] at 
t — > —00. Given the solution u((p,t), the amplitudes 
a„(t) can be found by the inverse Fourier transforma- 
tion: a n (t) = 5^ / Q 2lr u(ip,t)e~ lnip d(p. The solution of eq. 
(O) which obeys proper boundary conditions is: 



u((p, t) — exp 



-i 2g 




at \a\ 2 — 2/3, i.e. at 7 = 1 / ^ In 2 



3 ~ 0.254. Other two 
matrix elements (1/3/2)1/2,-1/2 = (^3/2)* 1/2,1/2 become 
zero at \a\ 2 = 1/3, i.e. at 7 = y^ln3 « 0.418. In the 
matrix U2 4 independent matrix elements have nodes: 
(U 2 )xx at \a\ 2 = 3/4; (U 2 ) 00 = at \a\ 2 = | (l±-*=); 

(U 2 )xo = at \a\ 2 = 1/2; (U 2 )x,-x = at \a\ 2 = 1/4. 

Electron motion driven by electric field. Here we 
present another exact solution. It describes transitions 
at intersection of infinite number of equidistant levels. 
This system can be described by the Hamiltonian, which 
matrix representation is 



tp ■ 



(18) 



Putting t = +00 in the solution (|18|) and taking the in- 
verse Fourier-transform, we arrive at the asymptotics: 



a n (t) « exp(— int /2 + i(n' — ri)ir/4) Ji„_ n /i (2V2%g) 



(19) 



Thus, the scattering amplitudes in terms of modified 
states with the fast phase factor exp(— int 2 /2) incorpo- 
rated are: 

(n I Uoo I ri) = e^ n '- n ^ 4 J ln _ n , { (2V2^g) (20) 
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It displays infinite number of oscillations with the LZ pa- 
rameter g. However, for large |n— n'\ the oscillations start 
with j > \n — n'\. For smaller values of g the amplitudes 
are small. 

This result can be easily extended to a more general 
non-diagonal matrix elements dep in which hopping from 
any site to any other site is allowed and its amplitude de- 
pends only on distance between sites. The corresponding 
Hamiltonian has the following matrix elements: 



fin 



9 k 



(21) 



For simplicity we present below the result for real g& 

g-k- 



(n\U 



i(n / — n)-7T / 4 

"hi 



2tt 



exp ( — i2v 2itf(ip, gj 



+i(n - n)ip) dip (22) 



where f(ip, gj) = V\. ^ cos kip. 

In conclusion, we presented a generalization of the LZ 
theory for arbitrary number of crossing equidistant lev- 
els in the case of time-dependent magnetic field and an 
electron in an infinite chain subject to a time-dependent 
electric field. In both cases high symmetry of the problem 
allows to find not only asymptotics, but also intermedi- 
ate values of the amplitudes. A number of analytical 
solutions exists for spin 1/2 in time dependent magnetic 
fields (l5) and all they can be generalized for higher spins. 
At the same time, our results are not universal. If the 
number of levels is more than two, they can approach 
to crossing point not equidistantly, and non-diagonal el- 
ements of Hamiltonian do not necessarily depend on one 
parameter only. In the majority of cases when the lev- 
els intersect systematically they remain equidistant, but 
the non-diagonal matrix elements of the Hamiltonian are 
different from above considered. Even in a simple case 
of a constant electric field perpendicular to the time de- 
pendent magnetic field, the Hamiltonian contains the ele- 
ments changing the projection m by ±2. For this case we 



have found general statement on separability and can re- 
duce the problem to already solved one for S < 2. On the 
other hand, our solutions demonstrate phenomena that 
are probably universal, like oscillations of the transition 
probabilities. 
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